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放射 治疗 〈 放 疗 ) 作为 一 种 治疗 恶性 肿瘤 的 治疗 手 
使 得 现代 放疗 又 迈 入 了 一 个 全 新 的 发 展 时 代 口 。 质 子 避 
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展 宽 布 拉 格 峰 (Spread-out Bragg Peak; SOBP) 
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质子 重 离子 治疗 计划 系统 利 
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描 仪 ， 相 对 阻止 本 领 


发 展 了 上 上 百年。 质子 与 重 离子 放疗 的 出 现 ， 
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HEHH CCT 作为 基础 输入 数据 ， 将 CT 值 (Hounsfield Unit, 
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与 SPR 或 WEPL 转换 曲线 。CT 图 像 以 每 个 像素 点 转换 成 WEPL 后 便 可 以 将 在 水 中 计算 得 到 的 剂量 映射 到 
CT 图 像 上 。 以 水 中 剂量 分 布 曲线 为 被 积 函数 以 每 个 像素 的 WEPL 为 积分 区 间 进 行 积分 再 除 以 相应 的 像素 
WEPL 即 可 得 到 每 个 像素 获得 的 剂量 53, 转换 曲线 的 准确 性 直接 影响 人 体 组 织 、 器 官 的 阻止 本 领 的 准确 性 ， 
从 而 影响 剂量 分 布 "。 获取 CT 值 相 对 阻止 本 领 转换 曲线 有 两 种 常用 方法 , 利用 动物 器 官 或 者 组 织 进 行 直接 
测量 "， 或 是 利用 组 织 等 效 模 体 进行 测量 得 到 *™。 上 述 方法 带 来 的 射程 不 确定 性 为 3%-5%， 需 要 在 治疗 计 
划 中 考虑 靶 区 外 扩 来 减 小 影响 ”。 


带电 粒子 穿 过 介质 时 ， 其 能 量 阻 止 本 领 可 用 Bethe-Bloch 公式 描述 ， 在 Schneider 等 中 的 研究 中 ， 介 质 


的 SPR 可 用 (1) 式 近 似 表达 ， 侈 去 了 公式 中 的 高 阶 无 穷 小 量 。 
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SPR ~ p (1) 


式 中 pe 为 介质 相对 于 水 的 电子 密度 ，B6 为 带电 粒子 的 相对 论 速 比 系数 ，mec? 为 电子 的 静 能 量 ，7 为 介质 材 
料 的 平均 激发 能 量 ， 水 的 平均 激发 能 量 [w=75eV™。 


从 (1) 式 可 以 看 出 SPR 除了 与 pe 呈正 比 关系 外 ， 还 与 介质 的 了 值 呈 一 定 的 相关 性 ， 而 SPR 对 6 的 依 
赖 性 非常 小 ， 低 于 0.4% 呈 ， 在 本 文 研究 中 忽略 此 项 影响 。 根 据 中 国 科学 院 近 代 物 理 研 究 所 研发 的 碳 离子 治 
疗 系统 使 用 能 量 范围 (120MeV/u 一 400MeV/) ,计算 时 统一 取 中 位 能 量 260MeV/u 时 的 8 值 ， 以 保证 最 


终 计 算 结 果 在 治疗 能 量 范 围 内 的 偏差 最 小 。 


在 Hünemohr 的 研究 中 ， 双 能 CT (DECT) 允许 提取 一 个 额外 的 参数 四， 从 而 解决 了 7 了 值 对 SPR 带 来 


的 不 确定 性 。 同 一 模 体 在 两 个 不 同 能 量 下 的 CT 值 允许 定义 除了 pe 之 外 的 每 个 体 素 的 有 效 原 子 序数 Ze, 


以 更 准确 的 信息 来 区 分 具有 相似 pe 的 组 织 ， 从 而 使 SPR 计算 更 为 精确 。 


本 研究 利用 CT 电子 密度 模 体 中 9 个 已 知 参数 的 组 织 等 效 模 体 插件 ， 根 据 Hinemohr 等 “在 研究 中 报 
道 的 DECT 双 参 数 提取 方法 ， 使 用 常规 单 能 CT (SECT) 设置 不 同 管 电压 对 组 织 等 效 模 体 进行 扫描 ， 提 取 
模 体 的 p。 与 Ze FIG SPR， 并 于 传统 方法 计算 的 SPR 做 对 比 ， 以 期 获 通过 此 方法 获得 更 加 精确 的 SPR 
值 。 


2 方法 与 材料 


2.1 设备 与 材料 


德国 西门 子 定位 CT 扫描 仪 ， 型号: SOMATOM Sensation Open, 美国 CIRS 公司 (www.cirsinc.com ) 


生产 的 电子 密度 模 体 ， 型 号 : 


L1 As 2s 


组 织 等 效 材料 


吸 气 肺 


呼 气 肺 


2.2 模 体 扫描 


EP 


Pe 


0.2 


0.496 


0.949 


0.976 


1.043 


将 组 织 等 效 模 体 插件 放置 在 配套 的 卡 槽 内 ， 使 用 定位 CT 进行 扫 


设置 为 80kVp、100kVp、120kVp 和 140kVp， 管 电流 统一 设置 为 320mAs， 


Z]v Ad 


厚 1.5mm. #4 


图 见 图 1. 


Zetr 


7.332 


7.382 


6.37 


6.766 


7.383 


CIRS 062M， 配 套 9 种 组 织 


等 效 模 体 插 


CIRS 062M 材料 参数 列表 


件 ， 参 数 如 表 1: 


组 织 等 效 材料 pe 
肝脏 1.052 
骨 松 质 200mg/cc 1.117 
骨 密 质 800mg/cc 1.456 
骨 密 质 1250mg/cc 1.695 


， 设 置 4 组 扫 


Zetr 


7.409 


10.55 


12.566 


13.456 


利用 RadiAnt DICOM Viewer 医学 影 


模 体 中 间 层 的 平均 CT 值 及 标准 偏差 ， 结 果 见 表 2. 4 组 序列 被 测量 


组 织 等 效 材料 
吸 气 肺 


呼 气 肺 


骨 密 质 1250mg/cc 


2.3 参数 提取 方法 


图 1. 模 体 扫描 示意 图 


像 浏览 器 软件 的 


扫描 层 厚 3mm, 


$i. 


表 2. 组 织 等 效 模 体 平 均 CT 值 /HU 及 标准 偏差 


80kVp 


-816.494 


-499.01+6. 


-97.53x5. 


-47.89+5. 


45.32+7. 


57.3847. 


355.36+ 


11.08 


1377.53£13.05 


2119.86+ 


16.26 


100kVp 


-817.69+ 


-500.0343. 


-85.12+3. 


-40.09+3. 


44.59+5. 


55.834 


4.55 


291.71+9.12 


1131.60+9.87 


1738.62+ 


11.21 


120kVp 
-818.12+5.81 
-500.33+3.28 
-76.68+3.06 
-34.3342.86 
45.16-44.26 
56.09+3.57 


255.78+7.52 


986.01+8.34 


1510.78+10.2 


140kVp 


-817.19+ 


-500.53+3. 


-72.72+2. 


-31.62+2. 


43.5724. 


53.9443. 


235.55+ 


888.444 


1343.12+ 


图 像 重 


"n 


站 序列 ， 管 电压 分 别 
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圆 形 选取 工具 ， 测 量 4 组 影像 序列 中 各 组 织 等 效 


又 域 空间 位 置 完全 一 


在 Hünemohr 的 研究 中 , 将 人 体 材料 的 X 冉 线 线 衰减 系数 Le 近似 分 解 为 一 个 虚拟 的 纯 光 电 吸 收 材 料 和 


一 个 虚拟 的 纯 康 普 顿 效应 材料 ; 


Zn 
Ug ~ Apefg + Bpezz (2) 


与 原子 序数 相关 系数 ， 左 是 关于 X 射线 能 量 近似 平坦 的 函数 ， 在 本 文 研 究 中 ， 


<a 


RP 为 光电 效应 截 
可 视 为 常数 。4、B 是 与 材料 几乎 无 关 的 比例 因子 ，Z 为 材料 原子 序数 ，E 为 X 射线 能 量 。(2) 式 对 CT 的 


X 射线 连续 谱 也 适 
对 于 CT 的 X 射线 连续 谱 ， 介 质 的 平均 线 误 减 系数 应 与 CT 图 像 的 FU 值 有 如 下 关系 : 


HU = (==) x 1000 (3) 
HE, 


其 中 ，FEw 是 水 的 平均 线 衰减 系数 。 后 续 计 算 中 ， 将 材料 相对 于 水 的 线 衰 减 系数 五 /KEw 记 为 4。 


通过 高 CHO R OO 两 档 能 量 扫 描 组 织 等 效 材料 ， 有 如 下 表达 式 : 


n = Apefg, + Bp . 
Hey ~ Apefg, + Bpe zur 
结合 (3) X WEBB A, Bla, K GO. (4) 式 ， 可 以 求解 出 p。 Zn， 并 将 已 知 数 据 代 入 ， 可 得 
出 如 下 结果 : 
pe = CX (I5 1) a -0 x (SR 1) 
Zerr = (em E 


FH (4) 式 可 知 ， 通 过 1 SO MAAS OA REE, aL RT DUE SC. D. XT Z-1-20 的 


范围 和 使 用 的 CT 光谱 ， 最 佳 系 数 为 n=3.1。 
在 本 文 研究 中 ， 使 用 全 部 9 种 材料 的 CT 影像 数据 ， 其 中 高 能 量 统一 采用 140kVp， 低 能 量 分 别 采用 


区 


80kVp、100kVp， 分 别 用 80/140 与 100/140 表示 不 同 组 别 ， 拟 合 出 C50 式 中 的 C、D 参数 ， 见 表 3。 


表 3. 参 数 C，D 拟 合 结果 


80/140 100/140 
C -0.8508 -1.6831 
D 5506.5149 10888.2124 


将 表 2 测量 结果 代入 (5) A, Hpo Ze HIR 1 的 已 知 参考 数据 做 比较 ， 结 果 见 表 4: 


d 4. 组 织 等 效 模 体 pe、2Qer 提 取 结 果 与 偏差 


DECT 


80/140 


-9.00% 
0.40% 
-0.11% 
0.61% 
-0.10% 
-0.10% 
1.52% 
1.10% 


-0.77% 


DECT 


100/140 


-8.00% 


0.60% 


-0.11% 


0.72% 


-0.10% 


-0.10% 


2.15% 


1.58% 


-1.06% 


UL X5 Zor 的 关系 如 下 : 


p. V TER Zr 计算 值 
DECT DECT DECT DECT 
组 织 等 效 材料 

80/140 100/140 80/140 100/140 
吸 气 肺 0.182+0.012 0.184+0.018 7.557+1.166 7.29542.438 
呼 气 肺 0.498+0.008 0.49940.011 7.535+0.376 7.506+0.497 
Aa 0.948+0.007 0.948+0.010 6.54540.229 6.49940.326 
乳腺 0.982::0.006 0.983+0.009 6.948+0.174 6.932+0.250 
肌肉 1.042+0.010 .042+0.014 7.497+0.197 7.504+0.312 
肝脏 1.051+0.009 .051+0.012 7.541+0.196 7.54940.265 
骨 松 质 200mg/cc 1.13+0.017 .141+0.025 9.665+0.195 9,52740.353 
骨 密 质 S00mg/cc 1.47240.022 .479:0.032 12.441+0.114 12.375+0.190 
骨 密 质 1250mg/cc 1.682+0.027 .677+0.039 13.57440.112 13.610+0.174 

2.4 In(D 值 计算 
C1) 式 中 的 In(D 值 使 用 上 述 方法 中 提取 的 Zer HEITI 

ee = 0.125Zeff + 3.379 Zeff < 8.5 


3 结果 


3.1 SPR 计算 


根据 表 4 的 计算 结果 ， 


结合 


CL» 


出 材料 的 相对 阻 上 本领， 并 于 参考 值 作 比较 ， 结 果 见 表 S. 


DECT 


80/140 


3.0796 


2.07% 


2.75% 


2.69% 


1.54% 


1.78% 


-8.39% 


-0.99% 


0.88% 


Zetr 偏差 


(6) 


中 SPR 参考 值 在 材料 的 说 明 中 并 未 直接 给 出 ， 表 中 的 值 是 表 1 中 提供 的 p。。Zr 值 ， 由 (1) 式 计 算 


表 5.SPR 测量 结果 与 误差 


组 织 等 效 材料 SPR 参考 

吸 气 肺 0.201 

呼 气 肺 0.497 

脂肪 0.965 

乳腺 0.987 

肌肉 .045 

肝脏 .054 

骨 松 质 200mg/cc .105 
骨 密 质 S00mg/cc .406 


SPR 测量 
DECT80/140 DECT100/140 
0.182+0.068 0.185+0.103 
0.498+0.017 0.499 +0.023 
0.961 +0.008 0.962 +0.011 
0.99 -- 0.007 0.991 +0.010 
1.042+0.010 1.042+0.014 
1.051 +0.009 1.051 +0.012 
1.134 土 0.015 1.142 土 0.022 
1.425+0.015 1.433 40.022 


DECT 80/140 


-9.49% 


0.19% 


-0.37% 


0.32% 


-0.24% 


-0.30% 


2.59% 


1.37% 


相对 误差 


DECT 100/140 


-8.15% 


0.43% 


-0.31% 


0.45% 


-0.25% 


-0.31% 


3.38% 


1.94% 


my 
Fu 


上 上 
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DECT 


100/140 


-0.50% 


1.68% 


2.03% 


2.45% 


1.64% 


1.89% 


-9.70% 


-1.52% 


1.14% 


骨 密 质 1250mg/cc 1.620 1.607 70.016 1.602 £0.023 -0.79% -1.13% 


32 5 SECT 单 参数 比较 


y = 0.0004x + 1.0183 
R’=0.998 HU>50 


y = 0.001x + 1.0177 
R? = 0.9977 HU<S0 
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2. 拟 合 HU-SPR 曲线 


采用 120kVp 扫描 参数 ， 使 用 传统 SECT 单 参数 提取 方法 "， 拟 合 出 HU-SPR 曲线 (图 2)， 再 通过 


线 ， 计 算 组 织 等 效 模 体 的 SPR 值 与 误差 ， 并 于 表 5 中 双 参 数 提取 方法 计算 SPR 的 误差 作对 比 〈 表 6、 


3)。 
表 6.DECT 方法 与 SECT 方法 SPR 测量 误差 结果 比较 
DECT SECT 
组 织 等 效 材料 

DECTS80/140 DECT100/140 SECT120 

吸 气 肺 -9.49% -8.15% -0.50% 

呼 气 肺 0.19% 0.43% 4.02% 
脂肪 -0.37% -0.31% -2.49% 

乳腺 0.32% 0.45% -0.41% 

肌肉 -0.24% -0.25% 1.72% 

肝脏 -0.30% -0.31% -1.23% 

骨 松 质 200mgy/cc 2.59% 3.38% 1.45% 
骨 密 质 800mg/cc 1.37% 1.94% 0.50% 


骨 密 质 1250mg/cc -0.79% -1.1396 0.1996 


4 讨论 


m SECT120 
* DECT100/140 


& DECT80/140 


p: 


- 10.00% 


- 5.00% 


ce EEE 


0. 


00% 


图 3. .DECT 方法 与 SECT 方法 SPR 测量 偏差 比较 


质子 重 离 子 放 疗 作 为 精准 放疗 手段 ,近年 来 , 引起 了 


转换 上 


线 对 于 质子 重 离子 放射 治 ; 


CT 设置 不 同 管 电压 时 产 


EET 


骨 密 质 1250mg/cc 


骨 密 质 800mg/cc 


质 200mg/cc 


国内 外 放疗 界 的 广泛 关注 。CT 值 与 相对 阻止 本 领 
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生 的 不 同 能 量 的 X 射线 对 组 织 等 


"m 


。 本 研究 通过 采用 常规 定位 


效 模 体 扫描 , 并 通过 高 低能 量 扫描 获得 的 两 组 CT 


值 ， 提 取 p。、&r 双 参 数 ， 并 通过 双 参 数 得 到 介质 的 相对 阻止 本 领 。 
本 文 方法 提取 的 参数 CR 4) pe 的 偏差 ， 吸 气 肺 在 DECT80/140 5 DECT100/140 组 的 计算 中 ， 分 别 


达到 了 -9.00% 与 -8.00%， 
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别 达 到 -8.39% 与 -9.70%， 其 余 偏差 均 小 于 396. Witt 等 ” 


据 进行 比较 时 ,发 现在 肺 部 等 效 模型 ， 
较 大 的 不 确定 性 所 。 因 


feet 


100/140 组 


孔 组 乡 
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间 上 绝对 


的 不 确定 
而 In(D 与 Ze BAER A, K 
出 显示 出 吸 气 肺 与 骨 松 质 的 1 
“ 均 偏差 0.12%。 其 结果 优 于 基于 SECT 单 组 
因此 ， 基 于 以 上 结果 ， 本 文 方法 能 更 加 精确 的 预测 SPR, 


的 相对 阻止 本 领 进行 进一步 校正 。 且 本 研究 使 
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匹配 ， 但 人 体 在 两 次 的 扫 


WE, Zrr 的 1 


此 ， 本 研究 中 ， 吸 气 肺 与 骨 松 质 


性 造成 的 。 在 Schneider[20 的 研究 ， 
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图 像 获 得 的 结果 C 
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由 于 器 官 运动 ， 会 带 来 
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! 均 在 190773. Zar 的 偏差 ， 骨 松 质 〈Bone200mg/cc) 在 两 组 的 误差 分 


在 研究 中 通过 理论 计算 得 到 的 数据 与 实际 测量 数 


差异 可 达 5%~15% , 其 分 析 主 要 是 因为 多 和 孔 材料 在 制作 过 程 中 
出 现 较 大 偏差 的 原因 ， 可 能 也 是 由 于 


孔 材 料 制 
，10% 的 7 值 误差 造成 的 SPR 预测 偏差 低 于 1.5%。 
来 的 SPR 预测 偏差 也 相对 较 小 。 表 5 的 SPR 计算 结果 中 ， 
i 差 较 大 ,除去 这 两 组 异常 的 数据 外 ,80/140 组 的 其 他 数 和 


的 可 行 性 ， 但 须 对 低 密 度 或 带 
j 的 是 固定 的 模 体 进行 两 次 CT H 
图 像 空间 的 不 匹配 。 
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EFH tit Zé 0.03%, 


， 其 两 次 图 像 在 
因此 ， 此 方法 只 能 


适用 在 一 些 短 期 内 〈 两 次 CT 扫描 期 间 ) 体位 重复 性 较 高 部 位 的 质子 重 离子 放疗 中 ， 比 如 使 用 高 精度 定位 
装置 固定 的 头 部 或 者 四 肢 部 位 。 其 他 部 位 如 使 用 此 方法 ， 只 能 使 用 DECT 扫描 ， 或 采用 形变 配 准 Pol， 使 两 
组 图 像 空 间 匹 配 。 

从 文中 (5) 式 分析 可 知 ， 不 同 能 量 组 的 计算 p。、Zr 差异 ， 首 先 由 高 能 量 的 HUj 值 决定 了 pe。、Zrr 的 基 
础 大 小 , 再 由 高 低能 量 的 CT 值 差 (HU — HUy) 在 前 者 基础 上 进行 调整 。 而 对 于 不 同 组 计算 中 , HU 的 值 都 


是 采用 的 最 高 能 量 140kVp 下 的 HU 值 ， 因 此, 决定 最 终结 果 精 度 的 主要 因素 是 高 低能 量 的 CT 值 差 。 基 于 


eu 


此 原因 ， 在 选择 此 方法 预测 SPR 时 ， 要 尽量 体现 出 CT 值 差 (EU -HUp)， 才 能 提高 预测 精度 。 因 此 ， 选 
择 能 量 差 比较 大 的 管 电压 设置 更 为 合适 。 

本 研究 的 SPR 参考 数据 ， 来 自 于 CIRS 062M 产品 说 明 ， 由 于 其 制作 过 程 中 的 不 确定 性 ， 本 身 精度 也 
均 有 一 定 的 偏差 。 因 此 ， 本 研究 获得 的 SPR 数据 需要 临床 实际 测量 进一步 校准 。 

此 外 ， 本 研究 中 拟 合 的 参数 数据 ， 仅 适用 于 实验 中 使 用 的 CT。 因 为 对 于 同一 种 组 织 等 效 模 体 ,使 用 不 
E CT 设备 、 不 同 扫描 参数 、 扫 描 时 模 体 所 在 位 置 以 及 模 体 周围 的 材料 等 均 会 对 CT 值 产 生 影 响 , 在 SECT 
方法 预测 的 SPR 中 ， 其 差异 可 以 达到 3%[21。 因此 ， 本 文 方法 也 有 必要 针对 不 同 的 CT 设备 ， 不 同 扫描 参 
数 制定 出 各 自 的 双 参 数 提取 方程 。 

总 之 ， 本 文 建立 了 设置 常规 CT 不 同 管 电压 获得 介质 的 双 CT 值 从 而 提高 对 介质 相对 阻止 本 领 精度 的 
方法 ， 对 进一步 提高 质子 重 离子 放疗 治疗 计划 剂量 计算 的 准确 性 具有 重要 参考 价值 。 


=r 
di 


` 


参考 文献 


[1] WANG Lan, DAI Xiaoya. Status and Prospect of Global Proton Heavy Ion Hospital [J].Hospital Construction and Equipment in China, 2016 
(1):27-31. Cin Chinese) 
CE, 戴 小 亚 . 全 球 质 子 重 离子 医院 现状 与 展望 中 . 中 国医 院 建筑 与 装备 , 2016(1): 27-31.) 


[2] CHUANG Y-C, WU P-H, SHEN Y-A, et al. Recent Advances in Metal-Based NanoEnhancers for Particle Therapy[J]. Nanomaterials, 2023, 
13(6): 1011. doi: 10.3390/nano13061011 
[3] LIU Xinguo, LI Qiang, DU Xiaogang, et al. Preliminary IMP Heavy Ion Therapy Planning System [J]. Nuclear Physics Review, 2010, 27 (4) : 
480-487. Cin Chinese? 
( 刘 新 国 , 李强 , 杜 晓 刚 , 等 . 初步 的 IMP 重 离子 治疗 计划 系统 [J]. 原子 核 物理 评论 , 2010, 27(4): 480-487.) 


[4 SCHAFFNER B, PEDRONI E. The precision of proton range calculations in proton radiotherapy treatment planning: Experimental verification of 
the relation between CT-HU and proton stopping Power[J]. Physics in Medicine and Biology, 1998, 43(6): 1579—1592. doi: 10.1088/0031-9155/43/6/016 
[5 LIU Xinguo, LI Qiang, DAI Zhongying. Dose calculation method in IMP heavy ion therapy [ J ].nuclear physics review, 2009,26 ( 1 ) : 69-75. 


Cin Chinese ) 


( 刘 新 国 , 李强 , 戴 中 颖 . IMP 重 离子 治 癌 中 的 剂量 计算 方法 [四 . 原子核 物 理 评 论 , 2009, 26(1): 69-75.) 


[6 MUSTAFA A A, JACKSON D F. The relation between X-ray CT numbers and charged particle stopping powers and its significance for 


radiotherapy treatment Planning[J]. Physics in Medicine and Biology, 1983, 28(2): 169—176. doi: 10.1088/0031-9155/28/2/006 


U RIETZEL E, SCHARDT D, HABERER T. Range accuracy in carbon ion treatment planning based on CT-calibration with real tissue Samples[J]. 


Radiation Oncology, 2007, 2(1): 14. doi: 10.1186/1748-717X-2-14 


8] SZYMANOWSKI H, OELFKE U. CT calibration for Two-dimensional scaling of proton pencil Beams[J]. Physics in Medicine and Biology, 
2003, 48(7): 861—874. doi: 10.1088/0031-9155/48/7/304 

9] COOLENS C, CHILDS P J. Calibration of CT Hounsfield units for radiotherapy treatment planning of patients with metallic hip prostheses: The 
use of the extended CT-Scale[J]. Physics in Medicine and Biology, 2003, 48(11): 1591—1603. doi: 10.1088/0031-9155/48/11/308 

0] YANG M, WOHLFAHRT P, SHEN C, et al. Dual- and multi-energy CT for particle stopping-power estimation: Current state, challenges and 
Potential[J]. Physics in Medicine & Biology, 2023, 68(4): OATROI. doi: 10.1088/1361-6560/acabfa 

1] SCHNEIDER W, BORTFELD T, SCHLEGEL W. Correlation between CT numbers and tissue parameters needed for Monte Carlo simulations of 
clinical dose Distributions[J]. Physics in Medicine and Biology, 2000, 45(2): 459—478. doi: 10.1088/0031-9155/45/2/314 

2] DEASY J. ICRU Report 49, Stopping Powers and Ranges for Protons and Alph Particles[J]. Medical Physics, 1994, 21(5): 709—710. Doi: 
10.1118/1.597176 

3] YANG M, ZHU X R, PARK P C, et al. Comprehensive analysis of proton range uncertainties related to patient Stopping-power-ratio estimation 
using the stoichiometric Calibration[J]. Physics in Medicine and Biology, 2012, 57(13): 4095-4115. doi: 10.1088/003 1-9155/57/13/4095 

4] MENG W-B, LI Q, YE Y-C, et al. Determination of virtual source position using back projecting zero and convergent arcTAN method for 
Scanning-passive scatter beam in carbon ion Therapy[J]. Journal of X-Ray Science and Technology, 2023, 31(1): 181—197. doi: 10.3233/XST-221274 

5] HÜNEMOHR N, KRAUSS B, TREMMEL C, et al. Experimental verification of ion stopping power prediction from dual energy CT data in 
tissue Surrogates[J]. Physics in Medicine and Biology, 2014, 59(1): 83—96. doi: 10.1088/0031-9155/59/1/83 

6] MURAISHI H, NISHIMURA K, ABE §, et al. Evaluation of Spatial Resolution for Heavy Ion CT System Based on the Measurement of Residual 


Range Distribution With HIMAC[J]. IEEE Transactions on Nuclear Science, 2009, 56(5): 2714-2721. doi: 10.1109/TNS.2009.2023520 


7] MENG Wanbin, MA Xiaoyun, ZHANG Yanshan, et al. Measurement of the CT-WEPL calibration curve in the carbon ion treatment planning 


system [ J ]. Modern Medicine and Health, 2022, 38 ( 6 ) : 920-924. (in Chinese) 


(Ji, SBS, KEL, 等 . 碳 离子 治疗 计划 系统 中 CT-WEPL 校正 曲线 的 测量 研究 [加 . 现代 医药 卫生 , 2022, 38(6): 920-924.) 


[18] WITT M, WEBER U, KELLNER D, et al. Optimization of the Stopping-power-ratio to Hounsfield-value calibration curve in proton and heavy 


ion Therapy[J]. Zeitschrift für Medizinische Physik, 2015, 25(3): 251—263. doi: 10.1016/j.zemedi.2014.11.001 

[19] URIE M, GOITEIN M, HOLLEY W R, et al. Degradation of the Bragg peak due to Inhomogeneities[J]. Physics in Medicine and Biology, 1986, 
31(1): 1-15. Doi: 10.1088/0031-9155/31/1/001 

[20] | LAO Zheng, FAN Wenhui, BI Fang, et al. Comparative study on deformation registration accuracy of medical image processing MIM software 
and Velocity software in elastic phantom [ J ].China Medical Equipment, 2021,18 ( 6 ) : 47-51. Cin Chinese) 

(Sill, BESCEE, HRI, 等 . 医学 图 像 处理 MIM 软件 和 Velocity 软件 在 弹性 模 体 中 形变 配 准 精度 的 比较 研究 [加 . 中 国医 学 装备 , 2021, 18(6): 47-51.) 


[21] SCHNEIDER U, PEDRONI E, LOMAX A. The calibration of CT Hounsfield units for radiotherapy treatment Planning[J]. Physics in Medicine 


and Biology, 1996, 41 (1): 111—124. doi: 10.1088/0031-9155/41/1/009 


Study on the method of acquiring the relative 
stopping power of charged particles based on dual- 


tube voltage CT images 
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Abstract : In proton and heavy ion radiotherapy, the relative stopping power of each voxel point obtained by CT 
scan for patient is the key to accurately calculate the dose distribution in the body using treatment planning system. 
In this study, a tissue equivalent phantom was scanned twice by setting different tube voltages for a conventional 
CT scanner. The relative electron density and effective atomic number of the phantom were calculated by the CT 
values derived from the two-set images, and then the relative stopping power of the phantom was calculated. The 
results showed that the relative stopping power of the tissue equivalent phantom acquired with the two-set CT values 
was better than that based on the conventional single set of CT images. Thus, the method established in this study 
might play an important role in further improving the dose calculation accuracy in the treatment planning of proton 
and heavy ion radiotherapy. 
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